Weak and Strong Localization in Low-Dimensional 
Semiconductor Structures 



S.-R. Eric Yang 1 ' 2 and J. Rammer 2,3 
1 Department of Physics, Korea University, Seoul 136-701, Korea 
2 Institute for Micro structural Sciences, National Research Council of Canada, Ottawa, Ontario, 

Canada, K1A OR6 

^Department of Theoretical Physics, Umed University, 901 87 Umed, Sweden 

(February 6, 2008) 

Abstract 



The dependence of the localization length on the number of occupied sub- 
bands N in low-dimensional semiconductors is investigated. The localization 
length is shown to be proportional to the number of occupied subbands in 
quasi-one-dimensional quantum wires while it grows exponentially with N in 
quasi-two-dimensional systems. Also a weak localization theory is developed 
for large N with a well-defined small expansion parameter 1/N. The temper- 
ature dependence of the conductivity deduced using this perturbation theory 

agrees with the experimentally observed dependence. 
PACS numbers: 71.55.Jv, 72.15Rn 
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According to the one-parameter scaling theory of localization electrons are localized 
in strictly one- and two-dimensional (ID and 2D) systems. The applicability of the weak 
localization theory of transport |2|] is thus limited in these systems: it is applicable when 
the Thouless length, Lt = \JDri n , is shorter than the localization length £ ||, and when 
a small expansion parameter exists (D and Tj n are the diffusion constant and the inelastic 
scattering time). Moreover, the expansion parameter of the diagrammatic transport theory 
{kpl) 1 ~ d H, where kp is the Fermi wave vector, is not small for d=l, and the localization 
length is comparable to the mean free path I. Recent advances in nanofabrication have made 
it possible to create quasi-one and -two dimensional (Q1D and Q2D) systems where several 
subbands may be occupied. The purpose of the present paper is to examine the range of the 
applicability of weak localization theory in Q1D and Q2D systems. The dependence of £ on 
the number of occupied subbands N determines the validity range of the condition Lt < £. 
When several subbands are occupied (fcpZ) 1_d is no longer the correct expansion parameter 
and a new small expansion parameter must be found. 

We apply Vollhardt and Wolfle's self-consistent diagrammatic theory of Anderson local- 
ization || to investigate these problems. Vollhardt and Wolfle's localization length for ID 
systems is in good agreement with the exact results of Berezinskii || and Abrikosov and 
Ryzhkin 0. In 2D their theory predicts that particles are always localized, in agreement 
with the one-parameter scaling theory of Abrahams et al. We extend their approach to 
the multisubband cases in the present paper. We find that in Q1D systems the localization 
length is proportional to N while in Q2D systems it grows exponentially with N. Our Q1D 
result is in good qualitative agreement with the previous results obtained using different 
methods ||. The basic physics behind the dependence of the localization length on N is 
the intersubband scattering which diminishes coherent backscattering leading to longer lo- 
calization lengths. We also find that 1/N{kpl) 1 ~ d serves as a small expansion parameter 
in Q1D and Q2D systems. Our Q1D result may provide an explanation for Mani and von 
Klitzing's experimentally found temperature dependence of the conductivity |||(T -p / 2 with 
p ~ 2). Their quantum wire has the number of occupied subbands between 7 and 11, and 
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therefore, according to our result, the localization length may be larger than the mean free 
path. So for temperatures where Lt is shorter than £ the conductivity may be evaluated 
using usual weak localization theory. 

Our interest is in conduction subbands of GaAs/AlGaAs quantum wells and wires and 
other subbands with similar electronic structures. In Q1D systems only the lowest subband 
originating from one of the confinement potentials is assumed to be relevant. Both in Q1D 
and Q2D systems all the relevant subbands are assumed to be parabolic with an effective 
mass m*, so each electron state is labeled by a subband index n and a wavenumber k. The 
level broadening due to the elastic scattering is assumed sufficiently small in comparison 
with the subband energy separations so that the subbands are well defined. In the presence 
of an external potential V ext (r) = V ext exp{iq - r — iut} with the wave vector q parallel the to 
x-axis, the Fourier transform of the induced density is given by <5n(q, ui) = x(q, u)V ext , where 
the Q2D and Q1D polarization functions are x(q, uj) = Jdzdz' x(q, z, z' ,uj), and x(q, uS) = 
fdydy'dzdz' x(<l,y,y' , z, z' respectively. The polarization function can be expressed in 
terms of the generalized diffusion coefficient D(q,u) [|TUJ 

X(q,w ) = N F — — — — 1 

— iuj + -L>(q, uj)q l 

where Np is the total d-dimensional density of states per spin at the Fermi energy. We use 
5-function impurity potentials Vi mp 5(r — r-), where r- are the positions of the impurities. In 
our model, where many subbands are occupied, the Fermi energy level broadening is assumed 
to be given by 7 = ttN^Uq with U$ equal to ni\Vi mp \ 2 (rij is the impurity concentration). 
We calculate corrections to the diffusion constant by including the weak localization ef- 



fects |n,I2|]. We include not only the intrasubband, but also the intersubband scattering. As 
long as the Thouless length is shorter than the localization length the correction to the con- 
ductivity due to weak localization effects is determined by summing the maximally crossed 
diagrams Green's functions are specified in the subband and momentum representation 
by respective set of quantum numbers. We have for the Cooperon 

A-m,n;m' ,n' s - + k , = ^m^m' ,n' 



+ J2 A^,„ ;M <IV (k + k', u) Ki, v , m ',n' (k + k', u) (2) 

i,v 

where 

n n ,„,(q, u) = X; G*(e + q - k)G#(e, k) (3) 

k 

Here G^(e, k) = (e — e n ^ + i/2r)~ 1 and the lifetime in each subband, r = l/2j, is approxi- 
mated to be independent of the subband index. The subband dispersion is e n ^ = e n + k 2 /2m 
where e n is the n'th subband energy, and wl v are the matrix elements of the impurity 
correlator Uq. 

The matrix equation Eq.(|2|) may be solved for A mjn;m / 5 „/ perturbatively in the small 
frequency, long wave length limit [13]. The matrix elements of the perturbation are 



V m ,n;m',n>{<\.,u) = A^.^, (IL n , n > (<l, iv) - 11^/(0,0)) (4) 

and the unperturbed eigenstates satisfy 

X A mim;M ,n M ,(0, 0)Sft, = a = 0, 1, 2, 3, ... (5) 

where 

n n ,n'( >°) = — t~ (6) 

e n ' - e n + i/T 

A perturbation calculation yields A m>m / ;ri)n /(q, uj) ,/AAo where the change in the 

eigenvalue Ao is given by 

AA =<$°|V|$ >=E<n'n„x(0,0) ( E Ky ;m y(k + k', w )^ im , (7) 

Inserting Eq.([|) and Eq.(||) into Eq.(^) we find 

a n,.^ ^_ -gT^Ayv; ^ -2^/irN F 

A„v4k + k,^- _ iw + Dc//(k + k , )2 ~ _^ + A) (k + k') 2 1 j 

where the effective diffusion constant -D e // is evaluated approximately as follows: Using 
^■mm-ii = j/nNp and replacing n nn (0, 0) with the average value < n nj „(0, 0) >= n < 
N n> F > /l ~ TfNp/N'y, we find the eigenvector [<3>° J 2 = j/ttNf- Again neglecting the 



off-diagonal elements of n n n /(0,0) we find an approximate value for the effective diffusion 
constant -D e // ~ D Q = J2 n D n /N m< v F > 2 r/d, where the prime indicate that we only sum 
over occupied subbands and < vp > is the average subband velocity on the Fermi surface. 
In our approximation each occupied subband has the same Cooperon and contribute to the 
conductivity independently. The weak localization correction leads then to a conductivity 
per spin of a Q1D system 

<n> Ne 2 r e 2 , r r „ I (L \. 

a = a + 5a = —(L — l) = cr 1 — — ( y — 1 ) 9 

The average density per subband <n> is < kp> /tt. From Eq.(|9]) we see that the correction 
due to weak localization effects is of order 1/N. When several subbands are occupied in 
Q1D and Q2D systems 1/N(kpi) d ~ l serves as the small expansion parameter. 

We now proceed to derive a formally exact expression for the inverse diffusion constant. 
We will then find an approximate expression for the inverse diffusion constant by replacing 
the irreducible vertex by the lowest order expression. The resulting equation is then extended 
into the insulator regime by demanding self-consistency. Our derivation is an extension of 
Vollhardt and Wolfle's method to QD systems. 

The impurity average of the product of Green's functions G^ k+ G^/ k / + may be expressed 

as 

^m,n,k;m' ) n' ) k'( C l; a; ) = ^mk + ^nk_^,k ,( 5n,n'Cm' + ^mk + ^rek_ rm,ji,k;m',n',k' ^^m'k' + ^Vk'- 

(10) 

The vertex function r minjk;m f in / k /(q, uS) is given in terms of the irreducible vertex function 

U m ,n,\i;m' ,n' , 

rm,nk;m>',k'(q)W) = fm,n,k;m',n' ) k'( c lj k>) 

+ X! ^m,n,k;m",n",k"(c[? ^)G^» k // + G^,, k ,,_ r m // >n // )k // ;m / )rl / >k / (q, u) (11) 
m",n",k" 

Using Eq.(O) $ may be written in terms of the irreducible vertex function 

$n,k;n',k'(<l> ^) = Cnk+^nk-^k'^n.n' + ^nk+^nk- X! ^n,k;n",k" (q ? <^)$n",k";n',k' (q, ^) (12) 

n",k" 



Using 

nR n A — AG n ^ . . 

Wk+^nk- - i ^ 6 

+ uj - k • q - AL„ ik 

we rewrite Eq.(p^) 

(uj - k • q - AE nik )$ n)k;n / )k /(q, u) = -AG n>k ^ k ,k"^n,n' 

+ f/ n,k;n",k"(q,^)$n",k»;n',k'(q,U;) (14) 

n",k" / 

Summing over n, k; n', k' in Eq.(14), and using AG nk = 2i^smG n ^ and the Ward identity 

As n>k = Y £4,k;<k'(q>^)A G <k' (15) 

n',k' 

we find 

u$(q,u) - q$j(q,uj) = 2niN F (16) 
where $(q,cj) = E n ,k;n',k' $n,k ; n',k'(q, w) and ^(q,u;) = E„,k ; n',k' ^?^,k,-»',k'(q, w). 



Multiplying Eq.([L4|) with k • q/k nF , where is the n'th subband Fermi wave vector, 



and summing over n, k; n', k' we find when 7 >> uj 

(q, ^)%(q, w) + i<?A)$(q, w) = (17) 

where 

K(iU,q) = l + T Yl ,2 Tt a -M -^G nik {U n>n .yy{k+k',Uj) - U ) AG n >y^ ■ q (18) 

n,k;n',k' ^Kp^npN 

In deriving Eq . ([17|) we have used the approximation AS njk = T, R k (E F + uj) — E^ k (E F ) ~ 
2i l 3mT l R k (E F ) = 2i~f, as << In addition we have used the peaked character of the 
spectral function to approximate 

$ n,k;ri',k'(q>^) ~ ^at^ AT I $ n',k';n",k" (q, w) 

n'.k' -2mN nF N y n , )k t^, )k „ 

dk-q^k'-q^^ \ 

+ -7 L 1? L $n',k';n"k»(q,^) (19) 

n',k' n",k" / 



keeping the first two terms in the Legendre expansion. Here we note that kq = k n pPi{cos 6), 
where 9 is the angle between k and q. 

Using the general form of the response kernel |J and approximating V n ,n;n,n — Uo by 

^n,n;n,n We find 

Do i+— y — - — -o (20) 

^ j\r_ ^ a, , i n r\2 v I 



L>(q, lu) ttN f ^ -iu + A)Q 2 

The residue of A ninjn) „ at iu = D^q 2 is 2 / ~f/irNF and is about 1/N times smaller than the value 
for the case where only one subband is occupied. This will have important consequences for 
the subband dependence of the localization length. Invoking self-consistency and multiplying 
Eq.([2DJ) with D(u)/D we find 

m = 1 _^_ T i (21) 

D nN F D q —iu/D(u) + Q 2 V ; 

This result determines the dependence of the localization length on the number of occu- 
pied subbands. By setting the upper Q-integration limit in Eq.(|2T|) to 1/7 and introducing 



a new variable y = JD(u)/ — iujQ we find 



"(") • ^' V ' d ; |- (22) 



D 7rN F D (27r) d ^ Jo 1+y 



where £ = y D(oo)/ — iool 2 = In the insulating regime the generalized diffusion constant 
approaches zero as u approaches zero ||15|| , and we find for Q1D and Q2D systems £ tan -1 £ = 
Ntt and £ = / (e N/x - l) 1/2 where A = 1/2ttE f t. For large N, we thus have £ = 2NI for 
Q1D systems. When N = 1 and Ep » r both results for Q1D and Q2D agree well with 
the perturbative estimates of £ in ID and 2D systems, namely irl and le nkp ^ 2 . We thus 
predict for quasi-dimensional systems a cross-over from the strongly localized to the weakly 
localized regime as the number of occupied subbands increases. The fact that the residue 
of A njn;n n is 1/N times smaller than that of singly occupied subbands is responsible for this 
effect. It is simple to show that the scaling functions f3 of Q1D and Q2D systems are not 
changed from the corresponding ID and 2D values. 

For sufficiently large N the condition < £ can be fulfilled and weak localization theory 
may be used to calculate the conductivity. We find in Q1D systems that 5cr ~ T~ l . The 



effect of Coulomb interactions in the presence of weak localization effects in ID and 2D 
systems has been investigated by Al'tshuler et al [JH]]. In Q1D and Q2D systems the same 
diagrams of Al'tshuler et al may be used when the Cooperon is replaced by our expression. 
The temperature dependence of these corrections is T~ 1 / 2 in Q1D systems. The temperature 
exponents found above, -1 and -1/2, are in agreement with the recent experimental values 
||. Mani and von Klitzing observed that the magnetoresistance of GaAs/AlGaAs quantum 
wires exhibits T _1 and T -1 / 2 behavior above 1.5 K, and at temperatures below 1.5 K the 
observed magnetoresistance saturates, as expected from the effect of electron localization. 
The number of occupied subbands in one of the wires used in the experiment is estimated 
to be between 7 and 11 fl7j . We believe that this large number for the occupied subbands 
is responsible for the observed temperature exponents. 
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